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We discuss the recent progress and the current status of experimental in- 
vestigations of spin-mediated energy transport in spin-chain and spin-ladder 
materials with antiferromagnetic coupling. We briefly outline the central 
results of theoretical studies on the subject but focus mainly on recent ex- 
perimental results that were obtained on materials which may be regarded 
as adequate physical realizations of the idealized theoretical model systems. 
Some open questions and unsettled issues are also addressed. 

PACS numbers: 75.40. Gb, 66.70. +f, 75.10.Pq, 75.^7.-m 

1. INTRODUCTION 

Heat transport in solids is often very sensitive to even weak disorder 
and is significantly influenced by phase transitions. That is why measure- 
ments of the thermal conductivity are an efficient experimental tool in solid 
state physics. Several books and review articles were devoted to thermal 
conduction in solids. '^'^'^ These previous authors usually discussed the heat 
transport by two species of itinerant quasiparticles only, i.e., by considering 
the thermal conductivity by phonons (^ph) and by itinerant charge carri- 
ers {tve)- Even when dealing with magnetic systems, the influence of mag- 
netic excitations is typically seen in a reduction of the phonon and electron 
thermal conductivities via scattering of these quasiparticles by magnons or 
magnetic impurities. The heat transport via the spin system itself was, to a 
large extent, neglected. One of the main reasons for this imbalance is that 
usually, it is very difficult to separate the spin contribution {ks) from the 
total measured thermal conductivity which, in most cases, is dominated by 
phonon heat transport. This problem is particularly severe for antiferromag- 
netic (AFM) three-dimensional (3D) materials where the linear dispersion 
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of the spin excitations causes, at low enough temperatures, Kg to adopt a 
dependence, the same as Atph(T). A much more favorable situation is met 
in strongly anisotropic spin systems, such as spin-chain materials, where 
the exchange interaction along a particular direction is much stronger than 
perpendicular to it. Since no thermal conductivity via spin excitations is 
expected perpendicular to the chain direction, investigating the anisotropy 
of the heat transport is an efficient method for the identification of Kg. What 
really makes spin-chain materials very attractive objects for transport mea- 
surements is, however, the influence of quantum effects on their properties, 
which are most pronounced for low spin values {S = 1/2 and S = 1). 

Unusual transport of energy and magnetization in certain ID quantum 
spin systems was predicted long agoP^I^ More recent theoretical studies 
along these lines induced a revival of interest in unusual transport proper- 
ties of integrable spin systems PI^'^I^ During the last decade, many new 
results have been obtained for transport properties of various ID quantum 
spin models, not only for the "simple" case involving an isotropic nearest- 
neighbor interaction, but also taking into account higher-order neighbor in- 
teractions, frustration effects and coupling to phonons. In addition, more 
complex ID spin arrangements such as ladders and zig-zag chains were con- 
sidered. The resulting theoretical progress was reviewed in recent articles 
I13|14[ This progress, linked to recent advances in preparation of materials 
whose features seem to fulfill some of the assumptions of the various ide- 
alized models quite well, has led to a considerable activity in experimental 
studies of heat conduction in ID spin systems. A corresponding summary 
of earlier activities was given in Ref. |T5l 

In this article, we present a brief survey of recent experimental results 
on the thermal conductivity via excitations in ID AFM spin systems. Heat 
transport in the 3D ordered spin state caused by weak yet unavoidable in- 
terchain interactions is not discussed here since it allows for a conventional 
description in terms of well-defined quasiparticles (magnons). We do not 
attempt to provide a full account of all the published papers on thermal con- 
ductivity of spin-chain materials, but rather concentrate on several central 
topics where some experimental progress has recently been achieved. First, 
we discuss whether the observed thermal transport in spin-chain compounds 
indeed reflects the expectations that emerged from the mentioned theoretical 
work. Then, we discuss what type of information on perturbations affect- 
ing the spin-mediated transport can be extracted from recent experimental 
data. Since most of the model materials are insulators, the dominating in- 
fluence is due to the interaction of the spin system with phonons as well 
as with magnetic defects. Another issue that we address here is how the 
predicted relations between thermal conductivity and spin conductivity in 
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low-dimensional spin systems were attempted to be verified in recent experi- 
ments. Finally, we discuss some features of the thermal conductivity of spin 
chains in external magnetic fields. 

2. THERMAL TRANSPORT IN SPIN CHAINS AND 

LADDERS 

2.1. Is it anomalous? 

In the modern theoretical literature, the energy transport in ID spin sys- 
tems has been addressed employing either the Boltzmann transport equation 
formalism or the linear response formalism. The former approach, which 
relies on the notion of quasiparticle modes with associated velocities and 
relaxation times, provides transparent results and is well suited for the anal- 
ysis of experimental data. However, this quasiparticle picture is not always 
applicable for quantum many-body systems. Thus, the most widely used 
approach was directed towards the calculation of transport coefficients via 
time-dependent current-current correlation functions, due to Kubo.'^ For 
the thermal conductivity at a finite frequency 

k{lo) = ^ r dte-'^' C'^ dT{3U-t - ir)jth), (1) 
J Jo Jo 

where jth is the energy current and (...) denotes the thermodynamic average. 
The real part of Eq. ([1]) can be decomposed into 

ReK(a;) = Dt\,5{uj) + Krcg(t^), (2) 

where the weight of the singular part (.Dth) is the so-called thermal Drude 
weight. The experimentally accessible quantity is the dc conductivity Re k,{ijJ - 
0). A nonzero -Dth implies a non-decaying energy current and, thus, bal- 
listic heat transport. Since the energy current is one of the conserved 
quantities ,'^'^1^ the existence of a nonzero Drude weight is often linked with 
the integrability of a system. The question whether Dth may, for various 
spin systems, adopt a nonzero value and how it evolves under the influence 
of different perturbations are the central points of current theoretical studies 
of thermal transport in quantum spin systems. The thermal Drude weight 
Dth{T) is as important for thermal transport as is the Drude weight D(T) for 
electric transport where its value at T = distinguishes an ideal conductor 
(-D(O) > 0) from an insulator(-D(0) = 0), as was put forward by Kohn.'^Sl Ex- 
tending Kohn's conjecture to heat transport and to nonzero temperatures,!^ 
several cases can be distinguished: 
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(a) thermal insulators if D^Y^iT) — and Kj.eg(w — > 0,T) = 0; 

(b) conventional thermal conductors if I?th(^) = and Kreg(i^ ~^ 0,r) > 0; 

(c) ideal thermal conductors if Dth (T) > for any value of 

However, an experimentally observed, anomalously large dc thermal 
conductivity does not necessary imply a ballistic thermal transport (case 
(c)). It is also possible that, even if DthiT) = (case (b)), the low-w 
part of Kreg(w,T) may adopt anomalously large values for certain systems 
with slowly decaying energy currents. Even if calculations for an idealized 
pure system suggest the existence of a nonzero Drude weight, in any real 
system, perturbations such as phonons, 3D couplings, or defects may destroy 
the integrability of the system. As a result, the 5-peak broadens into a 
Lorentzian with a width inversely proportional to the relaxation time t(T), 
the value of which depends on the type and strength of the perturbation, and 
the thermal conductivity is reduced accordingly. Thus, only a comparison 
of experimental data with theoretical calculations, taking into account both 
the singular and the regular parts of the conductivity spectrum, provides a 
reliable judgement of the situation. 

Real materials of current interest which may be regarded as reasonable 
physical realizations of model low-dimensional spin S systems belong to three 
classes: S = 1/2 chains, S = 1 chains, and two-leg S = 1/2 ladders. An 
often encountered case is the Heisenberg 5 = 1/2 XXZ chain, for which the 
Hamiltonian is 

H = J ^^{SfSf^i + SfS^j^^ + ASfSf^i), (3) 

i 

where J > is the intrachain nearest-neighbor exchange coupling and A 
characterizes the anisotropy of the interaction. The excitations are pairs of 
S = 1/2 quasiparticles (spinous) and the spectrum consists of a continuum 
with an upper and a lower boundary. It is gapless for |A| < 1 but a spin 
gap opens for larger |A|. The isotropic case A = 1 is especially impor- 
tant because a number of model systems with J ranging from quite high 
(e.g., SraCuOa with J/ks ~ 2000 K^pQl to rather low values (CuPzN with 
J/kB ~ 10 K^pS have been realized in appropriate material syntheses. Sys- 
tems which are described by Eq. ^ are integrable for |A| < 1 (easy plane) 
as well as for A = 1 (isotropic). In several theoretical papers and using dif- 
ferent methodical approaches it has unambiguously been demonstrated that 
at T > 0, the thermal Drude weight is nonzero if |A| < ij 22|23|24|2 5| other 
work established the thermal transport in S = 1/2 chains in the gapped 
phases, such as XXZ chains with Ising-type anisotropy |A| > Ijl^Slwith next- 
nearest neighbor interactions,'^ with three-spin interactions,'^ with dimer- 
ization and frustrations.'^S'SlESEo] T > 0, some calculations predicted a 
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nonzero Dth in the massive regime P^'^^'^ other studies question this con- 
clusion for dimerized and frustrated chains, howeverP^^ 

The same uncertainty about Dth{T > 0) exists for S = 1 chain systems. 
The excitation spectrum of the isotropic AFM S = 1 model has a large 
gap; the system is nonintegrable. Although earlier calculations provided 
support for a nonzero thermal Drude weight at all temperaturesp^ more 
recent work suggests that Dtii(T) vanishes in the thermodynamic limit.l^Slln 
the latter case, the thermal conductivity of the spin system is expected to 
be completely governed by the regular contribution Kj-cg{uj,T). 

A great deal of recent theoretical interest concentrated on the trans- 
port properties of two-leg S = 1/2 ladders. This interest was stimulated 
by experimental observations of a rather large spin thermal conductivity 
along the ladder direction in variants of (Ca,Sr)i4Cu2404i which contains 
spin ladders as a structural element.'^ The ladders are 5 = 1/2 chains (legs) 
connected in pairs, with an exchange constant J along the legs and J± per- 
pendicular to the legs (along the rungs). The interaction along the rungs 
is the integrability-breaking perturbation because for J± = 0, the system 
is just a set of 5 = 1/2 chains as described above. For J± > 0, the situ- 
ation with the thermal Drude weight is controversial, similar to other non- 
integrable systems with spin gaps mentioned above: the nonzero thermal 

Drude weight calculated for this system disagrees with more recent numeri- 
cal calculations. E l 3i | 34 | 35 | 

Measurements of the thermal conductivity were made for several phys- 
ical realizations of ID spin systems. Before discussing the experimental 
results for Kg, we wish to mention some aspects of ambiguity in the data 
analysis. The standard method of measuring the components of the thermal 
conductivity tensor k relies on the Fourier law of proportionality 

Jth = -kVT, (4) 

linking the heat flux Jth and the temperature gradient VT. The heat flux 
is typically directed along one of the main crystallographic axes, e.g. a, 
therefore we consider K^a = The measured thermal conductivity always 
includes contributions from all delocalized excitations acting as heat carriers. 
Therefore, separating the spin thermal conductivity Ks from contributions 
of phonons Kph, electronic quasiparticles Ke etc., is, in general, not straight- 
forward. Although Ke is negligibly small or zero for most studied spin-chain 
compounds, the phonon contribution is always larger or, at least, of the same 
order of magnitude as Kg. In principle, one may try to reduce Kph by intro- 
ducing some lattice defects which scatter phonons much stronger than spin 
excitations. This step has successfully been applied to the spin-ladder com- 
pounds (Sr,Ca,La)i4Cu2404i.'^ Unfortunately, the deliberately introduced 
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disorder produces some unavoidable additional scattering of spin excitations 
as well. 

The most widely used approach for establishing Kg exploits the fact 
that, because of very weak interchain interactions, the spin contribution to 
the transport perpendicular to the chains is usually negligible. Thus, for 
establishing the spin contribution one has to compare the thermal conduc- 
tivity measured parallel to the chains (which contains both Kph and Ks) 
with the thermal conductivity k"*" perpendicular to the chains (which con- 
tains only Kph). Still, because of the possible anisotropy of phonon transport, 
this method is, in some cases, only partly successful in avoiding substantial 
uncertainties in the evaluation of Ks{T). This is especially the case in the 
vicinity of phase transitions, where the phonon scattering by critical fluctu- 
ations is often very strong and not necessarily isotropic. For example, the 
existence of thermal conduction via the spin system was discussed for the 
S = 1/2 chain cuprate CuGeOa! ^^*^^ * ^^ * ^^ ^ Here the Cu atoms are linked via 
superexchange with Jc=10.4 meV along the chain direction (c axis), while 
in the perpendicular directions, the coupling is at least an order of magni- 
tude weaker (Ji, ~ O.IJ and Ja ~ — 0.01J).3IlThe temperature dependence of 
the thermal conductivity (T) along the chain direction exhibits a two-peak 
structure. It has been speculated that the higher-temperature peak is caused 
by the spin contribution to ^"(r), while the lower-temperature peak is of 
phononic origin .1^3351 However, the same features are also observed for k^[T), 
where Kg is expected to be negligible .'^S'^ Taking into account that the spin- 
Peierls transition occurs exactly between the two peaks at Tsp ~ 14 K, 
the two-peak feature may safely be attributed to phonon transport with 
enhanced scattering close to the transition.HSElESl As a result, despite the 
large amount of experimental work on thermal transport in CuGeOs, the 
identification of spin-carried heat conduction in this material has not been 
established unequivocally. A similar situation is met for the S = 1/2 chain 
compounds KCuFa and BaCu2Si207. Because of moderately weak interchain 
interactions, the 3D magnetic ordering sets in at relatively high temperatures 
(39.8 K and 9.2 K for KCuFa and BaCu2Si207, resp ectively) and is refiected 
in the temperature dependences of and k^.^SESI 

Only when the temperature dependences of the thermal conductivi- 
ties perpendicular to the chains are smooth and featureless, they can be 
used for evaluating the phonon background parallel to the chains. Exam- 
ples are shown in Fig. [TJ For the two materials, k^{T) exhibits a single 
peak at low temperatures and, at elevated temperature, k^{T) varies ap- 
proximately as which is typical for phonon thermal conduction of non- 
magnetic insulators.^ In turn, the thermal conductivity parallel to chains 
Ac'l(T) shows, on top of the phonon contribution, a distinct additional fea- 
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Fig. 1. (a) Thermal conductivity of the spin-ladder compound 

Sri4Cu2404i .'33133 The ladders are parallel to the c axis, (b) Thermal con- 
ductivity of the S = 1/2 chain material Sr2Cu03.4- The chains run parallel 
to the b axis. 

ture at higher temperatures. These extra contributions were attributed to 
thermal transport via the spin system, represented by Ks{T) = kI[T). The 
phonon background Kpj^(r) can be analyzed and subtracted from k"(T) but 
it is obvious that the uncertainties in Ks{T) are large for temperatures where 
the phonon background dominates. 

If we now turn to the experimental results so far obtained for Ks (T) , the 
main question is whether these results are consistent with the ballistic heat 
transport expected for integrable systems. If we define ballistic transport 
as transport where the mean free path of the excitations is limited by the 
sample dimensions, then the answer would be negative. The mean free path 
Is of ID spin excitations can roughly be evaluated using the kinetic equation 

Ks{T) = Cs{T)vs{T)h{T), (5) 

where Cs is the magnetic specific heat and Vg is the mean velocity of spin 
excitations. The data for ls{T) are available for S = 1/2 chain materials, 
such as Sr2Cu03,S3SrCu02,'32lBaCu2Si207,HslCuGe03,'™l39lYb4As3,IS21as 
well as for the S = 1 chain materials AgVP2S6,'^and Y2BaNi05,'^and some 
variants of the S = 1/2 two-leg ladder compound (Ca,Sr,La)i4Cu2404i.'32E6l 
In all cases, Ig was never larger than a few 10'^ A, much shorter than the 
lateral extensions of typical samples which were of the order of 1 mm. 

However, a more general definition of ballistic transport is provided 
by the presence of non-decaying contributions to the energy current which 
can be shown to be equivalent to the existence of a nonzero Drude weight. 
This definition is more appropriate because it does not rely on the valid- 
ity of the quasiparticle approach. Applying this criterion in analyzing the 
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experimental results of steady-state heat transport experiments, the identi- 
fication of the ballistic transport relies on demonstrating that the measured 
Ks exceeds the values of the thermal conductivity provided by the regular 
contribution in Eq. ([2]). Experimental data available in the literature indeed 
seem to support the theoretical predictions which claim ballistic thermal 
transport in 5 = 1/2 uniform Heisenberg chains and diffusive transport in 
S = 1 chains. In Fig. O we plot the spin-related energy diffusion constants 
De that were calculated from experimental data of Ks{T) for several spin- 
chain compounds^^ by using De{T) = Ks(T)/{Cs{T)a'^), where Cs{T) is the 
specific heat of the spin chain and a the lattice constant along the chains. 
For S = 1 chains, the energy diffusion constant has recently been calcu- 
lated from the regular contribution Kj-eg{u>,T) for high temperatures.'^ The 
high-temperature limiting value is indicated by an arrow in Fig. [2j For 
both S = 1 compounds presented in Fig. [2] the energy diffusion constants are 
clearly of the same order of magnitude as Dj^ values. In contrast, De{T) for 
the S = 1/2 chains are much larger than for the S = 1 species, they exhibit 
a strong temperature dependence and do not scale with T/ J, as one would 
expect for diffusion that is governed by intrinsic interactions only. Taking 
into account that the theory predicts Kj-cg{uJ,T) = for S = 1/2 chains in 
the gapless regime, this comparison strongly suggests that intrinsic diffusion 
is the dominating process in S" = 1 chains but much less so in S" = 1/2 
chains. A comparison between calculated values of Kreg{uj,T) and experi- 
mental results for Ks{T) has recently been made for two-leg spin ladders.!^ 
This comparison suggests that the experimentally observed thermal conduc- 
tivity via the spin system is governed by intrinsic diffusion. The emerging 
understanding of the dichotomy between diffusive and ballistic transport in 
ID spin systems is that the thermal transport in systems with a spin gap, 
such as S = 1 chains or S = 1/2 even-leg ladders, is determined by intrinsic 
diffusion, while for = 1/2 spin chains with gapless spectra, the thermal 
transport is ballistic and limited by external perturbations. 



Although most of the relevant theoretical work was based on the Kubo 
formalism, the analysis of experimental data on the thermal conductivity 
of low dimensional spin systems is usually employing Boltzmann's kinetic 
transport theory. In this approximation 
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Fig. 2. The energy diffusion constant De(T) calculated from the thermal 
conductivity data of the 5 = 1 chain compounds AgVP2S6 and Y2BaNi05, 
and of the 5 = 1/2 chain compounds BaCu2Si207 and SrCu02 (data from 
Refs. l46|49|51|52"]) . The arrow corresponds to the high-temperature limit 
calculated for 5 = 1 chains in Ref. [32l 

where /, e, Vg, and Is are, respectively, the distribution function, the energy, 
the velocity, and the mean free path of spin excitations with wavevector k. 
The mean free path Ig is related to the relaxation time t as Is = VsT. As- 
suming that several scattering mechanisms act independently of each other, 
in most cases a good approximation for the total mean free path is 

l;Hk,T) = J2l~l{k,T), (7) 

i 

where each ls^i{k,T) term corresponds to an independent scattering chan- 
nel. In order to keep the analysis tractable, in most experimental papers 
it is usually assumed that the scattering rates are /c-independent, such that 
ls,i{k,T) = ls^i{T). This assumption is often questionable but largely dic- 
tated by the lack of properly developed theoretical formalisms describing the 
relevant scattering processes of spin excitations in ID systems. Nevertheless, 
it is very useful in the sense that it allows to use the simplified expression 
of Eq. ([5]) to estimate the values and the temperature dependence of the 
average mean free path of excitations in the corresponding spin system. 

In many experimental publications on the spin thermal conductivity in 
ID spin chains and ladders,'^^ '^'^ ' ^'^ ' -'^" ' ^^ ' ^^ ' ^^ ' ^" ' ^^ ' ''^-'^ ' ^^ ' ^^! magnetic defects 
were considered as major scattering centers for spin excitations. It was 
always assumed that the corresponding scattering rates and hence the mean 
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free paths (/s,def ) were T- and /^-independent. In several studies, the influence 
of intentionally introduced disorder on the mean free path of spin excitations 
was investigated It was shown that the reduction of ls,dei indeed 
correlates with the concentration of magnetic defects, as is illustrated in 
Fig.Elfor the spin-ladder compound Sri4Cu2404i with Zn serving to partially 
replace Cu on the corresponding sublattice!^ Analogous correlations of the 
values of ls,def with the concentration of defects, estimated from NMR and 
magnetic susceptibility measurements, were reported in Refs. I48|52i 

There is no doubt that magnetic defects are very effective scatterers 
of spin excitations in ID magnetic systems. However, the convenient as- 
sumption that ls,dei is T- and fc-independent, seems not to be justified in 
general. The scattering of spin excitations by impurities and disorder was 
addressed theoretically in Refs. i58^59, for gapless S = 1/2 chains. The cal- 
culations suggest that the impurity relaxation time (and thus /s,def ) should 
be proportional to temperature. In Ref. [2^ the scattering by impurities 
was analyzed for ID spin systems with a gap in the excitation spectrum. It 
was found that the lifetime of spin excitations is energy dependent for the 
single-impurity case, and the situation is even more complicated for higher 
concentrations of impurities. 

Another species of scatterers needs to be taken into account in analyzing 
the data for the spin-ladder compound (Ca,Sr)i4Cu2404i where, depending 
on the chemical composition, different concentrations of mobile holes are 
present in the ladders. In Ref.[33lit was suggested that in Sri4_a;Cai.Cu2404i 
(x = 0, 2) the holes are the principal scatterers of magnons at high tem- 
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peratures. This conclusion was based on the observation that Is corre- 
lates with the average distance between holes. Subsequent experiments on 
(Ca,Sr)i4Cu2404i provided additional supporting evidence for the magnon- 
hole scattering in this materialP^'^ The holes are considered as a kind of 
mobile defects which scatter magnons in the same way as static defects, thus 
providing a T- and /c-independent term to the total inverse mean free path if 
the concentration of holes does not vary with temperature. At low tempera- 
tures, this assumption is not valid because the holes either adopt an ordered 
configuration or are transferred out of ladders.!^ 

At not very low temperatures, the scattering of spin excitations by 
phonons is expected to be significant. This scattering mechanism was taken 
into account in the analyses of experimental results for several spin-chain 
and spin-ladder materials by including a term ^7ph(^) right-hand 
side of Eq. ([7|) .ESESESBSES] Phenomenological expressions for ^7ph(-^) were 
suggested for spinon-phonon Umklapp processes in S = 1/2 chainpESland 
two magnon - one phonon scattering events in two-leg S = 1/2 ladders .'^Sl A 
theory for the spin thermal conductivity in S" = 1/2 chains, taking into ac- 
count spinon-phonon Umklapp scattering, was developed in Ref. [60l In this 
work, the transport in integrable systems with perturbations that weakly 
violate the conservation laws is considered. The presence of both spinon- 
spinon Umklapp processes and spinon-phonon scattering was demonstrated 
to lead to an exponential decrease of Hs{T) with increasing temperature, 
in apparent agreement with Ks{T) calculated from experimental data for 
Sr2Cu03. Spinon-phonon scattering in combination with impurity scatter- 
ing in S = 1/2 chains was theoretically analyzed in Refs. I58|59[ Three 
regimes with different temperature dependences for the thermal conductiv- 
ity were predicted: At low temperatures it is expected that Ks{T) oc T^. At 
intermediate temperatures, Ks{T) oc 1/T, and at high temperatures. Kg is 
expected to be constant. The experimental values for Ks(T) of Sr2Cu03, pre- 
sented in Ref. SHI seem to approach the predicted high-temperature constant 
value but, unfortunately, the temperature range covered by these measure- 
ments does not extend to high enough values in order to verify the predicted 
Kg = const satisfactorily. 

The role of phonons is not restricted to serving as a reservoir for the ex- 
change of energy and momentum with the quasiparticles of the spin system 
in the related scattering processes. First of all, a non-vanishing interaction 
between the lattice and the spin system is necessary for the very observ- 
ability of Kg in a typical thermal conductivity experiment, where the heat 
flux is first introduced into the phonon subsystem and then redistributed 
between phonons and spin excitations.'^ Since in an experiment probing the 
thermal conductivity, the spin excitations interact with phonons which are 
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not in equilibrium but have a nonzero net momentum along the temperature 
gradient, these interactions lead to a redistribution of momentum between 
phonons and spin excitations. If the scattering between phonons and spin 
excitations dominates over all other processes which lead to scattering within 
the spin systems (defects and Umklapp-processes), the resulting drag effects 
may modify Kg ■ The influence of such drag effects on the thermal conductiv- 
ity of conventional (3D) FM and AFM systems was discussed in Refs. I62|63[ 
Depending on the relative strengths of magnon-phonon, magnon-defect, and 
phonon-defect scattering, and on the type of defects (magnetic or nonmag- 
netic), complex k{T) variations were predicted. Recently, the problem of 
spin drag has theoretically been addressed for spin ladders .El It was shown 
that spin-phonon drag effects modify both Kph and Hs and also cause interfer- 
ence effects which are taken into account by the introduction of Hs,ph- With 
respect to experiment, no comparison of experimental data with spin-phonon 
drag predictions have been published up to now. 



2.3. Influence of a magnetic field 

According to recent theoretical work, an external magnetic field can 
modify Hs not only because it changes the spectrum of magnetic excitations 
but also via so-called magnetothermal corrections.'^^^ In general, the heat 
current Je and the magnetization current Jm are related to the temperature 
gradient VT and the magnetic field gradient V-B by 

\ _ ( Lmm Lme \ ( Vi? I , s 

3E )~\ Lem Lee ) \ -VT ) ' ^''> 

In zero magnetic field, Lem = Lme = and hence, the coefficient Lee '^^ 
equivalent to the thermal conductivity Kg (in the following only the magnetic 
contribution to the heat transport is considered). In magnetic fields, the off- 
diagonal elements in ([8]) may be nonzero, and this leads to a number of in- 
teresting effects. For example, by analogy with the Seebeck effect in electric 
conductors, where the heat flux induces an electric voltage along the direc- 
tion of the flux, jE is expected to induce a gradient of magnetic induction in 
magnetic systems. A number of interesting predictions for magnetothermal 
effects in 5 = 1/2 chain compounds which fulfill the XXZ-model conditions, 
can be found in recent work in Refs. l65|25|67[ At present, no experimental re- 
sults on magnetothermal effects in spin chains are available, mainly because 
of considerable practical difficulties in carrying out such experiments.'^ 

Nonzero off-diagonal elements in ([8]) also modify the thermal conduc- 
tivity Kg of the spin system .'^SEsMl Since heat transport experiments are 
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normally done under the condition of Jm = 0, it follows from the definition 
of K in Eq. (g]) that, for B ^ 0, 



The second term on the right-hand side of Eq. ([9]) represents the mag- 
netothermal correction, which was theoretically investigated in detail in 
Refs. I25|661 Nontrivial dependences of these corrections on temperature, 
magnetic field and the anisotropy parameter A in ([3]) were predicted. Un- 
fortunately, an experimental evidence for magnetothermal corrections is still 
missing. Most of the good physical realizations of the ID spin model sys- 
tems discussed above, such as Sr2Cu03, are not suitable for this type of 
experiments, because of the extremely large exchange constants with J/kB 
of the order of 10^-10^ K. Available laboratory magnetic fields are much too 
small to influence such systems in any significant way. Moreover, Shimshoni 
et al^ argued against the existence of these corrections in real materials 
because the conservation of the total magnetization may be broken by, e.g., 
spin-orbit coupling. In their paper, also the influence of magnetic field on 
the spin thermal conductivity ol S = 1/2 chains was addressed. Complex 
fractal-like dependences of Ks on the magnetization, resulting from field- 
induced variations of spin-phonon Umklapp scattering, were predicted. 

The authors of Ref. 66 studied the relation between the thermal Drude 
weight Dth and the spin Drude weight Kg 



which is analogous to the celebrated Wiedemann-Franz law for electrical con- 
ductors. The behavior of the parameter Lq in different parts of the phase 
diagram of the = 1/2 XXZ chain was investigated. The relation between 
the heat conductivity and the spin conductivity was tested experimentally 
in studies of the 5 = 1/2 chain compound Sr2Cu03 and the 5 = 1 chain 
compound AgVP2S6, respectively. In Fig. |H we show the comparison of the 
energy diffusion constant De{T) obtained from the zero-field thermal con- 
ductivity measurements ,'3SE] and the spin diffusion constant Ds{T) estab- 
lished by NMR experiments.'SS'^ For Sr2Cu03, both transport parameters 
show similar temperature dependences and differ by a factor close to 2. Also 
for AgVP2S6, De and Ds adopt similar absolute values and exhibit very 
similarly shaped temperature dependences. Earlier theoretical worliP^*^ for 
5=1/2 and 5=1 chains and the recent calculations for 5 = 1 chains,'351 
indeed, predict values of the ratio De/Ds in the interval between 1.4 and 3 
at high temperatures. 



f^s = Lee 
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Fig. 4. The energy diffusion constants De{T) (open symbols) of Sr2Cu03 
and AgVP2S6, estimated from thermal conductivity dataj^SEU and the spin 
diffusion constants Ds{T) (solid symbols) from NMR measurement^^^*^ 
(from Ref.lSI]). 



3. SUMMARY AND OUTLOOK 

Recent experimental studies accumulated a considerable amount of ev- 
idence for heat transport via magnetic excitations in paramagnetic one- 
dimensional spin arrays. In some cases, particularly for 5 = 1/2 spin chains, 
the spin thermal conductivity is rather large, consistent with theoretical 
predictions which claim a ballistic form of transport. Magnetic defects were 
identified as one of the important sources for the scattering of spin exci- 
tations. The experiments also suggest that spin-phonon scattering is very 
effective in reducing the energy flow, especially at high temperatures. Re- 
cent measurements provided some support for diffusive heat transport in 
ID spin systems exhibiting a gap in the excitation spectrum. Nevertheless, 
the situation is far from being clear, analogous to a similar controversy in 
the theoretical sector. Special analyses of experimental data revealed a link 
between the heat conduction and the spin transport in spin-chain systems. 

In spite of the experimental progress in certain number of prob- 

lems could not yet be solved. One of the central issues in current theoretical 
studies is whether the transport in spin chains with a gap in the excitation 
spectrum is ballistic or diffusive. The available experimental data sets do not 
provide enough solid evidence in favour of either of these expectations. The 
number of suitable materials that were studied so far, is simply too small. 
For example, only two S = 1 chain compounds were investigated. In both 



Thermal conductivity in spin-chain materials 

materials, the concentration of magnetic defects was relatively high, leaving 
the possibility that intrinsic features of the thermal conductivity are masked. 
Also for the probed spin-ladder materials, the experimentally observed large 
heat flow carried by spins does not seem to contradict the possibility of dif- 
fusive transport governed by intrinsic interactions. In order to clarify this 
situation, more experimental work is needed. 

Concerning the integrable spin systems with gapless excitation spectra, 
the most important problem is to establish, how the spin thermal conductiv- 
ity is affected by various possible integrability-breaking perturbations, such 
as dimerization, frustration, weak interchain coupling, next-nearest-neighbor 
interactions, anisotropic interactions, various types of disorder, and coupling 
to phonons. Very little experimental work has been done in this direction. 
Even for the most extensively studied case of scattering by defects, the T - 
and /c-dcpcndence of the corresponding scattering rate is not yet established 
unambiguously. For progress in this direction, model materials that allow 
for a controlled alteration of different perturbations are required. 

In spite of the rich variety of phenomena predicted by theory, very lit- 
tle experimental work has been done on the influence of magnetic fields on 
the spin-carried heat conduction in ID spin systems. Progress in this area 
requires materials with weak intrachain exchange and yet much weaker inter- 
chain interactions. Besides variations in the thermal conductivity, a number 
of interesting new features in magnetothermal transport are expected to be 
observed. 
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